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VIBRATION  TEST  OF  RICHARD  B.  RLSSELL  CONCRETE 
DAM  BEFORE  RESERVOIR  IMPOUNDMENT 

PART  I:  INTRODUCTION 

1.  Information  from  forced  vibration  tests  of  Richard  B.  Russell  Dam 
before  impoundment  of  the  reservoir  was  used  to  determine  the  natural 
frequencies,  mode  shapes,  modal  damping  ratios,  and  relative  joint  motions. 
This  low-level  forced  vibration  test  conducted  in  January  and  February  1982 
seized  a  rare  opportunity  to  experimentally  measure  the  dynamic  properties 
of  a  concrete  gravity  dam  without  hydrodynamic  interaction.  The  experi¬ 
mentally  measured  dynamic  properties  can  be  compared  to  2-D  and  3-D  finite 
element  dynamic  modeling  to  check  the  validity  of  these  analyses.  A  forced 
vibration  test  planned  in  1984  after  impoundment  of  the  reservoir  will  pro¬ 
vide  an  experimental  measure  of  the  prototype  hydrodynamic  interaction  and  a 
comparison  of  the  changes  of  the  dynamic  properties. 

2.  Richard  B.  Russell  Dam  has  recently  been  built  by  the  Corps  of 
Engineers  approximately  170  miles  (274  km)  from  the  mouth  of  the  Savannah 
River  between  the  States  of  Georgia  and  South  Carolina.  A  view  of  the  dam 
from  the  downstream  Georgia  side  is  shown  in  Figure  1.  As  shown  in 


Figure  1.  View  of  Richard  B.  Russell  Dam  from 
downstream  Georgia  side 
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Figures  2  and  3,  the  crest  of  the  concrete  gravity  part  of  Richard  B.  Russell 
Dam  is  1884  ft  (574  m)  long  and  is  composed  of  13  nonoverflow,  8  intake,  and 
11  spillway  monoliths,  the  tallest  of  which  is  approximately  200  ft  (61  m) 
high. 

3.  Another  role  of  measuring  the  dynamic  properties  of  the  dam  is  to 
provide  a  record  of  the  dam's  structural  condition  as  built.  If  the  dam's 
structural  integrity  is  suspect,  following  damage  sustained  from  an  earthquake 
or  other  sources,  the  results  of  test  data  before  and  after  the  damage  would 
be  useful  in  determining  the  significance  of  the  damage  (Richardson  1980) . 
Finally,  the  results  from  the  first  vibration  test  also  can  be  used  to  evaluate 
dynamic  analyses  previously  conducted  to  assure  the  safe  seismic  design  of 
the  project  (Norman  and  Stone  1979;  Norman  1079). 
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Figure  2.  Plan  and  elevation  views  of  Richard  B.  Russell  Dam 


PART  II:  EXPERIMENTAL  PROCEDURE 


4.  Richard  B.  Russell  Dam  was  excited  by  the  device  shown  in  Figure  4. 

It  consists  of  a  17,000-lb  (7,711-kg)  inertial  mass,  which  was  driven  by  a 
Zonic  ZTL  electrohydraulic  serve-controlled  actuator.  The  shaker  was  welded 
to  a  4-ft  by  5-ft  by  2-in. -thick  (1.22  by  1.52  by  0.0508  m)  steel  plate  which 
was  epoxied  to  the  mass  concrete  of  the  dam.  The  epoxy  used,  w  Jika's 
Sikastix  370  Sikadur  HiMod  and  it  performed  very  well  at  the  c  temperatures 
encountered  during  the  curing  process.  Three  plates  were  affi.  at  the 
exciter  locations.  The  force  input  to  the  dam  was  computed  as  roduct  of 

the  measured  acceleration  and  the  17,000-lb  mass.  A  typical  fo  ..  spectrum 

of  the  shaker  for  the  tests  is  shown  in  Figure  5.  The  maximum  stroke  of  the 
actuator  is  limited  to  1  in.  (0.0254  m) .  From  25  to  6  Hz  (1  Hz  =  1  Hertz 
=  1  cycle/second),  the  force  spectrum  is  flat  (for  a  peak  force  output  of 
20  kips  (89  kN)).  The  rolloff  from  6  to  1  Hz  is  due  to  the  stroke  limitation. 

5.  Q-Flex  Model  No.  QA1000  servo-accelerometers,  as  shown  in  Figure  6, 
with  a  sensitivity  of  0.25  volts/g  were  used  to  measure  the  horizontal  ac¬ 
celerations  of  the  concrete  dam.  The  accelerometers  have  a  frequency  response 


Figure  4.  Inertial  mass  vibrator  -  located  on 
monolith  7 


range  from  0  to  500  Hz.  Analog 
signals  were  recorded  by  a  Sangamo 
Sabre  V  tape  recorder  on  32-track 
magnetic  tapes.  The  tape  recorder 
has  a  frequency  response  range 
from  0  to  2,500  Hz.  A  Structural 
Dynamics  Analyzer  HP5423A  manu¬ 
factured  by  the  Hewlett-Packard 
Company  was  used  to  later  digitize 
and  process  the  analog  signals. 

6.  All  the  recording 
instrumentation  was  housed  in  a 
closed  van  positioned  on  the 
crest  of  the  monolith  adjacent  to 
the  drive  point  as  shown  in 
Figure  7a.  The  generator  and 
the  hydraulic  power  supply  for 
the  shaker  are  seen  to  be  lo¬ 
cated  on  the  same  monolith  with 
the  van.  The  contractor's 
gantry  crane  was  required  to  place 

all  the  experimental  equipment.  . 

Figure  6.  Accelerometer  mounted  on 

In  Figure  7b  the  van  is  seen  monolith  22 

being  lifted  to  the  crest  of  the  dam. 

All  the  equipment  is  listed  in  Table  1.  A  schematic  of  the  excitation  setup 
is  shown  in  Figure  8. 

Vibration  Tests 

7.  The  three  exciter  locations  were  chosen  at  monoliths  7,  16,  and  22. 
These  positions  roughly  correlate  to  the  sixth,  half,  and  quarter  points  of 
the  dam,  respectively.  Monolith  7  is  the  tallest  nonoverflow  section  and  the 
existence  of  vertical  shafts  in  monoliths  16  and  22  made  it  convenient  for 
placing  instrumentation  for  monitoring  the  cross  section's  horizontal 
response.  In  Figures  3  and  9  the  cross  sections  and  the  locations  of  the 
drive  point  monoliths  are  shown.  The  force  input  in  the  downstream-upstream 


and  setup  at  monolith  22 


b.  Lifting  of  van  to  crest 


direction  to  the  dam  was  computed  as  the  product  of  the  measured  acceleration 
and  the  17,000-lb  mass  (7,711  kg).  The  three  exciter  locations  were  at 
monoliths  7  (elevation  494),  16  (elevation  483),  and  22  (elevation  494). 

8.  Ideally  it  would  be  desirable  to  excite  a  structure  at  one  location 
and  mount  enough  transducers  at  all  points  of  interest  to  measure  the  response 
of  the  structure.  However,  on  a  structure  as  large  as  Richard  B.  Russell  Dam 
this  was  not  possible  because  of  the  number  of  available  transducers  and  the 
construction  on  portions  of  the  dam.  The  first  series  of  tests  to  measure 
horizontal  accelerations  of  the  dam  in  plan  and  cross  section  began  with 
drive  point  7  (monolith  7).  The  first  array  of  accelerometers  was  placed 
along  the  crest  at  the  center  line  of  monoliths  1  to  11  to  monitor  the 
horizontal  upstream-downstream  motion.  The  accelerometer  locations  are  shown 
in  Figure  9  for  all  the  crest  tests.  The  second  and  third  arrays  of  crest 
accelerometers  were  placed,  respectively,  at  the  center  line  of  monoliths  8  to 
16  and  16  to  22  for  drive  point  16.  The  final  two  arrays  of  crest  accelerom¬ 
eters  were  placed  at  the  center  line  of  monoliths  16-22  (complete  overlap  of 
the  third  crest  array)  and  22  to  31  for  drive  point  22. 

9.  At  each  drive  point  the  corresponding  monoliths'  accelerometers 
were  placed  to  measure  the  distribution  of  horizontal  acceleration  with  eleva¬ 
tion.  The  location  of  accelerometers  with  respect  to  elevation  of  the  mono¬ 
liths'  cross  sections  is  shown  in  Figure  10. 

10.  Relative  joint  motion  arrays  were  positioned  at  three  locations  and 
consisted  of  two  accelerometers  closely  spaced  on  each  side  of  the  monoliths' 
joint  being  monitored.  The  joint  arrays  were  at  elevation  483  joint  15-16 
and  elevation  436  spillway  crest  joints  21-22  and  22-23.  The  joint  acclerom- 
eter  locations  are  shown  in  Figure  11. 

11.  There  were  eight  distinctive  test  arrays  and  sixteen  separate  tests 
were  conducted.  All  the  tests  were  1/2*  decade/minute  sine  sweeps  from  25  Hz 
to  1  Hz.  This  rate  was  found  to  be  slow  enough  for  a  required  quasi-steady 
state  condition  to  provide  enough  data  in  a  given  time  interval  for  adequate 
digital  time  resolution.  Each  test  lasted  about  2.8  minutes;  but  there  were 
many  hours  of  preparation  required  to  conduct  a  test  for  a  given  array.  The 
electrical  signals  output  from  the  accelerometers  and  the  shaker  were  recorded 
as  analog  signals  on  the  32-track  tape  recorder. 


e  (refe 
ites . 


sweep  from  20  Hz  to  2  Hz  (1  decade) 


Figure  10.  Cross  sectional  accelerometer  locations 


IEX9I 


Figure  11.  Relative  joint  accelerometer  locations 


Dynamic  Cylinder  Tests 


12.  Dynamic  tests  were  conducted  on  5  of  22  standard  b-  by  12-in.  (152 
by  305-mm)  cylinders  cast  from  concrete  mixtures  used  in  the  construction  of 
the  dam.  The  cylinders  were  placed  in  a  fog  room  within  48  hours  of  casting. 
Before  being  transported  to  WES  the  cylinders  were  removed  from  the  fog  room. 
Once  at  WES  the  cylinders  remained  in  the  testing  room  up  to  test  time. 

13.  For  the  dynamic  tests  the  cylinders  were  supported  horizontally  at 
the  node  points  (corresponding  to  the  fundamental  flexural  mode  of  vibration 
of  a  free-free  beam,  see  Figure  12).  With  an  accelerometer  mounted  on  a  test 


CABLE  TO 
ANALYZER 


CONCRETE 
6-BY  12-IN. 
TEST  CYLINDER 


FORCE-INSTRUMENTED  HAMMER 


—  DIRECTION  OF  MEASURED 
ACCELERATION 


DIRECTION  OF 
MEASURED  FORCE 


AXIS  OF  CYLINDER 


FUNDAMENTAL 
MODE  SHAPE, 

OF  CYLINDER 
ALONG  ITS  AXIS 


SOFT  SUPPORTS  (AT  NODES 
OF  THE  FUNDAMENTAL 
MODE  SHAPE) 


Figure  12.  Schematic  of  dynamic  cylinder  test  setup 


cylinder  a  force-instrumented  hammer  was  used  to  strike  the  cylinder  with 
the  applied  force  parallel  to  the  measured  acceleration.  This  method  of 
testing  has  provided  reasonable  results  in  a  number  of  similar  applica¬ 
tions  (Alexander,  A.  Michel,  1981). 
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PART  III:  DATA  REDUCTION  AND  EXPERIMENTAL  RESULTS 


14.  To  obtain  the  necessary  data  form  for  the  modal  analysis  the  HP5423A 
Structural  Dynamics  Analyzer  was  used.  The  HP5423A  has  a  two-channel  digital 
processor  with  12-bit  resolution  per  channel.  All  the  analog  signals  were 
played  back  by  using  the  force  input  and  one  of  the  accelerometer's  output  for 
that  test  to  digitize  an  input/output  (I/O)  pair.  A  total  of  70  1/0  pairs 
were  digitized  representing  the  transfer  functions  (or  frequency  response 
functions)  at  the  58  accelerometer  locations  on  the  dam.  A  transfer  function 
is  the  complex  (i.e.,  composed  of  real  and  imaginary  parts)  ratio  of  the 
measured  output  acceleration  to  the  measured  input  force.  The  data  were 
stored  on  HP  data  cartridges  compatible  with  the  HP5423A.  Besides  the 
frequency  response  functions  (or  transfer  functions)  the  coherence  functions 
from  strategic  points  were  also  computed  and  saved.  The  coherence  function 

for  an  I/O  pair  measures  what  ratio  of  the  output  acceleration  is  linearly 
caused  by  a  measured  input  force:  thus  each  I/O  pair  represents  a  single 
input/single  output  model.  For  a  more  detailed  description  of  spectral 
functions  an  excellent  reference  is  Engineering  Applications  of  Correlation 
and  Spectral  Analysis  (Bendat  and  Piersol  1980) . 

15.  Transfer  and  coherence  functions  are  shown  in  Appendix  A  which 

are  representative  of  the  overall  response  of  the  dam.  The  coherence  function 
is  nondimens ional  and  has  values  from  0  to  1  (0  to  100%  coherent  signal). 

The  natural  frequencies  were  determined  by  a  quadrature  peak  picking  method 
(called  a  line  cursor  fit  in  the  HP  users  guide  (Hewlett-Packard  Co.  1981  and 
Hewlett-Packard  Co.  1982)).  This  involves  observing  the  peaks  in  the 
imaginary  part  and  the  zero  crossings  in  the  real  part  of  the  transfer 
function  for  a  particular  frequency.  The  HP5423A  fits  a  curve  through  11 
data  points  of  the  peak  in  the  imaginary  plot  of  the  transfer  function  and 
computes  the  modal  damping  ratios  and  shape  parameters.  Therefore,  one 
transfer  function  is  representative  of  the  overall  structural  damping  for  a 
particular  resonant  frequency  and  mode  shape.  Ten  mode  shapes  were  in¬ 
vestigated  and  the  modal  parameters  are  tabulated  in  Table  4.  For  an  explana¬ 
tion  and  example  of  the  modal  parameter  extraction  method  see  Appendix  B. 

16.  Before  the  modal  analysis  was  performed  a  geometry  of  the  dam 
represented  by  the  accelerometer  locations  was  described  and  the  coordinates 
were  entered  into  structural  component  tables  in  the  HP5423A.  Using  the 


estimates  of  natural  frequencies  and  damping  ratios  residue  values  are  com¬ 
puted  from  each  transfer  function  of  the  response  points.  The  residues  are 
used  to  compute  shape  vectors  for  each  response  point.  Then  the  shape  vectors 
combined  from  each  response  point  define  the  mode  shapes  of  the  dam  for  each 
resonant  frequency  investigated.  The  10  identified  mode  shapes  are  shown  in 
crest  plan,  three  cross  sections,  and  isometric  views  in  Figures  13-22. 

17.  The  relative  joint  motions  are  determined  by  faking  the  ratio  of 
one  transfer  function  to  the  transfer  function  across  the  joint.  Thus  the 
input  term  cancels  out  resulting  in  the  ratio  of  accelerations  across  the 
joint.  The  relative  joint  motions  are  shown  for  joints  15-16,  21-22,  and 
22-23  in  Figures  23-28. 

Dynamic  Cylinder  Tests 


18.  The  dynamic  cylinder  test  results  were  used  to  independently 
measure  material  damping  and  the  dynamic  modulus  of  the  concrete  used  in  con¬ 
structing  the  dam.  The  values  of  damping  obtained  from  the  dynamic  cylinder 
tests  are  a  lower  bound  that  could  be  expected  in  the  dam  structural  system. 
Using  the  HP5423A  the  electrical  signals  from  the  accelerometer  and  the  load 
cell  (force  instrumentation)  were  processed  and  a  transfer  function  was  com¬ 
puted.  From  the  transfer  function  the  fundamental  resonant  frequency  was 
determined,  and  using  an  appropriate  frequency  bandwidth  a  damping  ratio  was 
computed.  The  results  are  listed  in  Table  3. 
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Figure  16.  Mode  shape  4  (Sheet  1  of  3) 
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Figure  18.  Mode  shape  6  (Sheet 


view 


Figure  19.  Mode  shape  7  (Sheet 
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Figure  22.  Mode  shape  10  (Sheet  1  of  3) 
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Figure  24.  Coher 


Figure  26.  Magnitude  of  the  transfer  functions  across  joint  22-23 
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Figure  28.  Ratio  of  accelerat 


PART  IV:  DISCUSSION  OF  RESULTS 


19.  In  interpreting  the  foregoing  results  it  should  be  recognized  that 
the  conditions  of  the  dam  deviated  slightly  from  the  experimentally  optimal 
situation.  Earth-moving  equipment  was  being  operated  24  hours  a  day,  and  there 
was  water  flowing  through  the  sluices  constantly  even  though  there  was  no  ap¬ 
preciable  reservoir  water  impounded.  Finally,  there  was  on-going  construction 
changing  minutely  the  characteristics  of  certain  sections  of  the  dam  between 
tests.  However,  the  data  were  adequate  to  conduct  a  thorough  modal  analysis. 

In  Table  2  the  water  levels  (headwater  and  tailwater)  are  listed  for  certain 
tests.  These  three  activities  probably  contributed  some  random  excitation  to 
the  dam  in  addition  to  input  from  disturbances  of  nonstructural  components 
caused  by  the  controlled  input  to  the  dam  during  the  tests.  The  abutment  on 
the  Georgia  side  was  not  completed  which  makes  the  west  edge  of  the  crest  of 
monolith  1  a  "free  end."  The  existing  grade  of  the  earthen  dam  on  Georgia's 
abutment  monolith  1  was  approximately  60  ft  (18.3  m)  below  the  crest  eleva¬ 
tion  495.00.  South  Carolina's  abutment  was  finished,  restraining  monolith  32 
with  an  expected  degree  of  fixity. 

20.  The  three  mode  shapes  at  5.94,  6.75,  and  7.56  Hz  in  Figures  13,  14, 
and  15  correspond  to  the  first  three  natural  modes  of  vibration  found  through 
test  and  analysis  of  the  Pine  Flat  Dam  (Rea,  Liaw,  and  Chopra  1972).  Some 
discontinuities  appear  in  these  lower  mode  shapes.  In  Figure  5  notice  that  the 
force  rolls  off  at  these  frequencies.  With  a  lower  force  input  a  lower 
signal-to-noise  ratio  (SNR)  results.  A  low  SNR  implies  the  desired  signal 

is  more  contaminated  by  noise  than  a  measurement  with  a  higher  SNR.  The  SNR 
can  be  calculated  from  the  coherence  function  for  any  frequency  as 

SNR  =  -1-  l - 

1  -  Y  (f) 

2 

where  y  (f)  is  the  amplitude  value  of  the  coherence  function  at  frequency  f  . 
For  instance  in  Figure  13  the  SNR  at  10  Hz  is  only  (0.60)/(l  -  0.60)  =  1.5:1. 
Thus  the  apparent  anomalies  in  the  shapes  could  be  attributed  to  noisy  data. 

An  undetected  natural  mode  of  vibration  may  exist  between  7.6  and  8.1  Hz 
since  mode  No.  4  in  Figure  16  identified  at  8.19  Hz  resembles  in  crest  plan  a 
fifth  normal  mode  shape  of  vibration  of  a  fixed-fixed  beam  (Shock  and  Vibration 
Handbook,  ed .  Harris  and  Crede,  1961).  It  is  possible  that  this  mode  was  not 
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picked  up  because  the  three  vibrator  locations  were  nodes  of  this  vibrational 
shape,  or  very  close  mode  coupling  occurred  with  the  third  or  fifth  mode,  such 
that  the  resonance  for  the  fourth  mode  is  overshadowed.  This  mode  coupling 
(very  closely  spaced  resonant  frequencies)  is  prevalent  above  8  Hz.  The  mode 
shapes  at  9.06,  10.06,  and  10.94  Hz  seem  to  correspond  to  the  sixth,  eighth, 
and  tenth  normal  mode  shapes  in  crest  plan  as  seen  in  Figures  17,  18,  and  19 
of  a  fixed-fixed  beam.  Above  11  Hz  the  spatial  density  of  response  points  is 
too  sparse  to  narratively  characterize  the  natural  mode  shape  for  the  entire 
dam.  The  response  of  the  dam  above  11  Hz  appears  to  be  characterized  by 
resonances  of  the  different  monolith  sections.  For  instance  the  spillway  sec¬ 
tions  have  large  responses  at  13.12,  14.62,  and  16.37  Hz  but  the  intake  sec¬ 
tions  are  not  as  responsive  at  these  frequencies,  as  seen  in  Figures  20,  21, 
and  22.  Similarly  the  intake  and  nonoverflow  sections  are  responsive  at  some 
frequencies  when  the  rest  of  the  structure  is  not.  Some  nonlinear  responses 
are  attributed  to  the  monoliths'  joint  behavior.  The  analysis  of  the  relative 
joint  motions  reveals  this  nonlinear  behavior.  From  observing  the  crest  plan 
view  of  the  mode  shape  at  5.94  Hz  it  appears  that  the  spillway  sections  are 
the  most  stiff  and  the  intake  sections  the  most  flexible.  By  observing  the 
first  three  mode  shapes  (in  Figures  13,  14,  and  15)  the  estimated  location  of 
the  center  of  stiffness  is  between  monoliths  11  and  13.  The  structural  damp¬ 
ing  ratios  for  the  modes  shown  range  from  1.43  to  5.16  percent.  These  values 
bound  the  damping  ratios  determined  from  the  Pine  Flat  vibration  tests  (Rea, 
Liaw,  and  Chopra  1972). 

21.  For  most  of  the  frequencies  the  governing  mode  shape  in  cross 
section  resembles  the  fundamental  mode  shape  of  a  cantilever  beam  (Part  b 
of  Figures  14-22).  Only  for  monolith  16  does  an  intermediate  node  exist 
above  13.13,  as  shown  in  Figures  20b,  21b,  and  22b.  Consistent  with  recent 
conclusions  about  dam-foundation  interaction  (Chopra,  Chakrabarti,  and  Gupta 
1980)  there  is  motion  detectable  near  the  base  of  each  monolith  and  as  ex¬ 
pected  is  less  than  the  response  at  the  crest.  The  measured  mode  shapes  in 
the  cross  sectional  views  are  also  consistent  with  those  from  the  Pine  Flat 
vibration  test  (Rea,  Liaw,  and  Chopra  1972). 

22.  The  instrumented  joints  16-15,  22-21,  and  22-23  were  monitored 
close  to  the  input  since  these  joints  would  have  the  most  relative  motion.  At 
joint  16-15  in  Figures  23  and  24  the  transfer  function  and  the  coherence 
function  help  explain  what  physically  occurs.  The  peaks  (resonant  frequencies) 
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correlate  well  in  both  transfer  functions,  and  above  10  d  tin. 
function  is  mostly  above  0.90.  The  comparison  of  the  transfer  functions  re¬ 
veals  that  there  is  some  relative  motion  at  resonant  trev'  iu  i«.-s  but  ;  he  slip 
is  very  small.  The  larger  relative  motion  indicated  by  •  r.-.i  ■.  ot  ’he  ac¬ 
celerations  of  joint  16  to  15  (A,, /A,,)  in  Figure  5  o-u.  •  •  1  ;  •  urs  it  non¬ 
in  lb 

resonant  frequencies.  The  joint  may  be  acting  as  a  tilt.  ■  di  ,-,i-  'in  the 
energy  at  these  frequencies.  At  frequencies  where  A  / A  ; -  .  i ■  •  itet  than 
1.0  also  defines  the  ratio  of  two  very  small  mot  ions  rec  t  •  tin  resonant 

responses.  Where  the  coherence  function  drops  snoot  !  ■-  •  d  ■  . , ?  •  •.  c  in¬ 
dicate  nonlinear  behavior  of  the  joint  or  the  preset'..  •  ■  .•  r-.  SNK  as 

seen  below  10  Hz.  The  value  of  A,./Alr  equal  to  i  .(•  i  ":p .  . .  •  •,  ■  'el.it  ive 
joint  motion. 

23.  The  joints  22-21  and  22-23  should  be  expectec  ' o  behave  similarly 
as  these  joints  are  both  at  a  spillway  crest.  In  Figures  2"  and  17  the 
transfer  and  coherence  functions  are,  respectively,  shown  tor  joint  22-23. 
Again  the  peaks  occur  at  the  resonant  frequencies.  Below  10  Hz  the  coherence 
functions  are  decreasing  indicating  decreasing  SNR.  However,  at  17.0  and  17.8 
Hz  there  are  two  smooth  drops  in  the  coherence  function  for  joint  side  23, 
and  the  transfer  functions  reveal  large  differences  in  responses.  This  may 
indicate  constant  relative  joint  motion  above  15  Hz.  The  two  nulls  in  the 
coherence  function  at  17.0  and  17.8  Hz  correlate  with  the  indication  of  two 
frequencies  where  joint  slippage  is  most  excited.  The  ratio  of  A^^/A^^  shows 
a  similar  plot  in  Figure  28  where  above  15  Hz  there  is  the  greatest  measure  of 
relative  joint  motion. 


Discussion  of  Dynamic  Cylinder  Tests 

24.  The  dynamic  testing  of  the  standard  6-  by  12-in.  (152-  by  305-mm) 
concrete  cylinders  provided  a  measure  of  the  material  damping  of  the  dam's 
concrete.  In  Table  3  it  is  seen  that  the  material  damping  ratios  range  from 
0.315  to  0.537  percent.  The  damping  ratios  of  the  dam  determined  from  the 
modal  analysis  represent  the  combined  effects  of  boundary  and  material  damping 
Thus,  the  effects  of  the  boundary  conditions  (the  supports:  foundation  and 
abutments)  are  greater  than  the  effects  of  material  damping.  This  is  con¬ 
sistent  with  a  previous  analytical  study  (MacBain  1972). 


Comparison  to  Previous  Analysis 


25.  Previously  a  2-D  dynamic  finite  element  (FE)  analysis  was  conducted 
on  the  modified  geometry  of  the  nonoverflow  monolith  to  compute  the  funda¬ 
mental  frequency  (Norman  1°79).  The  analysis  assumed  no  reservoir  with  the 
height  of  the  monolith  at  185  feet.  A  value  of  7.63  Hz  was  obtained  com¬ 
pared  to  5.94  Hz  identified  as  the  fundamental  frequency  of  the  prototype 
structure.  This  difference  is  expected  and  the  value  of  the  2-D  FE  analysis 
would  differ  for  the  following  reasons:  (1)  Foundation  interaction  is  not 
accounted  for  in  the  2-D  FE  analysis.  The  foundation  interaction  decreases 
the  frequency  (Norman  and  Stone  1979).  (2)  The  elastic  modulus  in  the  2-D  FE 
analysis  was  a  static  value.  Under  low  level  dynamic  conditions  the  modulus 
increases,  thus  both  the  stiffness  and  the  frequency  increase.  (3)  The  3-D 
FE  behavior  is  shown  to  be  stiffer  than  the  2-D  FE  behavior  of  Pine  Flat  Dam 
for  the  fundamental  frequency  and  the  frequency  increased  (Rea,  Liaw, and  Chopra 
1972).  The  3-dimensional  response  is  expected  to  differ  from  the  2-dimensional 
finite  element  analysis.  In  the  case  of  Pine  Flat  Dam,  the  3-D  response  was 
stiffer  in  the  fundamental  mode  than  the  2-D.  The  effect  in  Russell  Dam  is 
not  necessarily  the  same  since  Pine  Flat  and  Russell  Dams  do  not  have  geo¬ 
metrically  similar  boundary  conditions. 

26.  Other  factors  would  also  cause  a  difference.  One  factor  also  not 
accounted  for  in  the  2-D  model  is  the  relative  joint  motion  between  monoliths. 
Relative  joint  motion  is  a  nonlinear  behavior  and  it  may  increase  the  flexi¬ 
bility  of  the  system. 

Structural  Integrity  of  Richard  B.  Russell  Dam 


27.  As  discussed  previously,  the  dam's  natural  frequencies,  damping 
ratios,  and  the  mode  shapes  are  consistent  with  design  analyses  and  the 
dynamic  properties  determined  from  the  Pine  Flat  prototype  vibration  test. 
Therefore,  the  dam  appears  to  be  structurally  sound  as  built. 


PART  V:  CONCLUSIONS 


28.  The  low  level  vibration  test  of  Richard  B.  Russell  Dam  has  shown 
that  the  electrohydraulic  shaker  is  adequate  to  excite  the  structure  with 
measurable  vibrations.  Also,  the  use  of  a  structural  dynamics  analyzer  is 
a  very  efficient  tool  to  both  digitize  and  analyze  the  data  to  compute  the 
dynamic  properties. 

29.  The  natural  frequencies,  damping  ratios,  and  mode  shapes  determined 
are  reasonably  consistent  with  design  analyses.  The  natural  frequencies  for 
the  10  mode  shapes  investigated  expectantly  ranged  from  5.97  to  16.4  Hz.  The 
damping  ratios  ranged  from  1.43  to  5.16  percent  which  bounds  the  values  from 
the  Pine  Flat  vibration  tests  (Rea,  Liaw,  and  Chopra  1972). 

30.  The  first  seven  shapes  of  the  dam  resemble  in  crest  plan  the  first 
three,  fifth,  sixth,  eighth,  and  tenth  normal  mode  shapes  of  a  fixed-fixed 
beam.  The  governing  mode  shape  for  monoliths  7  (tallest  nonoverflow),  16 
(spillway  transition),  and  22  (spillway)  at  the  resonant  frequencies  in¬ 
vestigated  resembled  the  fundamental  mode  shape  of  a  cantilever  beam.  The 
tallest  section,  monolith  16,  also  revealed  a  mode  shape  at  14.6  and  16.4  Hz 
resembling  the  second  normal  mode  of  a  cantilever  beam. 

31.  Relative  joint  motions  exist  and  are  a  function  of  the  excitation 
frequency.  Similar  behavior  was  observed  between  joints  16-15  and  the  spill¬ 
way  crest  joints  22-21  and  22-23.  There  appears  to  be  a  filtering  effect  be¬ 
tween  resonant  frequencies  at  joint  16-15  and  little  relative  movement  at  the 
resonant  frequencies  below  15  Hz.  However,  above  15  Hz  (actually  15-20  Hz) 
there  is  larger  relative  motion  between  monoliths  21,  22,  and  23  at  the  spill¬ 
way  crest  than  at  joint  16-15.  The  nonlinear  behavior  of  the  dam  is  attributed 
mostly  to  the  monoliths'  joint  behavior. 

32.  Since  the  dynamic  properties  (natural  frequencies,  damping  ratios, 
and  mode  shapes)  determined  are  reasonable  and  contain  no  gross  anomalies,  the 
dam  appears  to  be  structurally  sound  as  built. 
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PART  VI:  FUTURE  WORK 
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33.  Dynamic  finite  element  analyses  are  planned  to  study  the  2-D  and  3-D 
behavior  of  the  dam.  The  "as  tested"  prototype  structural  conditions  are  to  be 
modeled.  These  dynamic  analyses  will  check  their  validity  when  compared  to  the 
prototype  test  results.  The  foundation  interaction  effects  studied  will  be 
used  also  as  a  check  and  compared  to  the  prototype  test  results. 

34.  The  mode  shapes  that  have  been  obtained  will  provide  good  compari¬ 
sons  to  2-D  and  3-D  and  dynamic  finite  element  modeling. 
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35.  The  second  vibration  test  of  the  dam  will  be  conducted  with  a  full 
reservoir.  From  Table  2  it  seems  also  the  temperatures  during  tests  were 
moderate  which  implies  there  is  no  need  to  conduct  comparison  vibration  tests 
only  during  the  month  of  February.  However,  when  a  vibration  test  is  con¬ 
ducted  with  the  reservoir  impounded  it  is  more  desirable  to  choose  a  period 
when  the  pool  water  elevations  will  not  change  significantly.  At  that  time 
it  may  be  possible  to  determine  mode  shapes  of  certain  sections  by  just  measur¬ 
ing  the  accelerations  caused  by  the  ambient  noise. 
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36.  The  second  test  will  be  analyzed  with  a  dynamic  FE  model  including 
the  hydrodynamic  interaction.  Also,  the  second  vibration  test  will  provide 
another  measure  of  the  dam's  dynamic  properties  after  the  reservoir  is  filled. 
Hence  largely  the  effect  of  hydrodynamic  interaction  will  be  measured. 
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Table  1 

Equipment  List 


Accelerometer 


Amplifiers 


Tape  recorder 
record 
playback 


Electrohydraulic  vibrator 
(horizontal  mode) 


Hydraulic  power  supply 


Generator 


Digital  processor 
and  modal  analyzer 


SDC ,  Q-FLEX  Model  QA1000 

sensitivity,  vol ts/g . . . . 0 . 25 

Built  by  ISD ,  WES 
phase  compensated 


Sangamo  Sabre  V 

Sangamo 

Zonic,  ZTL  electro  hydraulic 
servo  control 

Driving  Mass . . . . 17 ,000  lb-mass 
(7,711  kg) 

Dynamic  force  (pk  sinusoidal)  lb- 
force  used  -  20,000  lb  (88.96  kN) 

Frequency  range,  HZ  1  -  25  HZ 
total  harmonic  distortion  25% 

Operating  pressure  3000  psi 
(20.67  MPa) 

Pump  capacity  -  70  gpm  driven  by  a 
150  hp,  3  phase,  440  v  motor 

Diesel  Caterpillar  3306T,  150  KW  set, 
Ser.  No.  66D41617 

HP  5423A  Structural  Dynamics 
Analyzer 
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Recorded 

Temperatures 

and  Water  Levels 

Du  ring  Tes t  s 

Temperature 

Test 
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Wa  1 1'  r  h  1  ev.i  i 

i  Ui ,  t  t 

- 

No. 

Date 

Test  Time 

Monol i th 

Headwa t  e i 

'.i  1  1  V*  .1 1  ( ‘  i 

1 

1/30/82 

67 

7 

__ 

.. 

j 

2 

1/30/82 

60 

7 

-- 

-- 

3 

2/1/82 

50 

7 

-- 

-- 

■ 

5 

2/7/82 

50 

7 

349 . 8 

3  iO  .  8 

1  ■ 

6 

2/7/82 

52 

7 

349.8 

330.8 

m 

7 

2/10/82 

43 

16 

342.0 

332.0 

‘ 

1 

8 

2/13/82 

60 

16 

340.9 

331.9 

» 

9 

2/15/82 

60 

16 

340.9 

332.0 

■ 

1 1 

2/19/82 

55 

22 

352.5 

332.5 

- 

12 

2/20/82 

69 

22 

_  _ 

_  _ 

% 

15 

71 

22 

333.0 

332.7 

2/21/82 

16 

2/21/82 

71 

22 

333.0 

332.7 

Average  = 

342.7 

331.9 

U 

±7.48 

±0.76 

1 

i 

Coefficient  of 

Cy  =  0.02 

0.002 

variation : 

1 

'  • 

« 

l  • 

« 

.  j 

Table  3 

Frequency  and  Damping  Measurements  on  6"  by  12"  Concrete 
Cylinders  from  Richard  B.  Russell  Dam 


Des 

ign 

Damp  i  rig 

Strength 

Cylinder 

Frequency 

Dynamii 

Mo  tin  1  u  s 

Ra  t  io 

ksi 

KP 

No. 

Position-' 

HZ 

x  1 0°  jj  s  l 

x  10°  KPa 

0/ 

/o 

4.0 

27.6 

5037 

1 

4023 

5 . 59 

38.5 

0.331 

1 . 0 

20.7 

5464 

1 

3953 

5 . 40 

37.2 

0.3  30 

5464 

2 

3978 

5.46 

37.6 

0.315 

3.0 

20.7 

5166 

1 

4167 

5.82 

40.  1 

0.537 

5166 

2 

4133 

5.84 

40.2 

0.391 

5166 

3 

4088 

5.77 

39.8 

0.339 

2.0 

13.8 

4954 

1 

3643 

4.58 

31.6 

0.366 

4954 

2 

3611 

4.50 

31.0 

0.380 

2.0 

13.8 

5833 

1 

3625 

4.54 

31.3 

0.364 

Cylinder 

No. 


Cylinder 

No. 


Results  from  Uniaxial  Compression  Test 
i  ri...  •  Uniaxial 


Poisson' s 
Ratio 


Initial  Elastic 
Concrete  Modulus 


10  psi 
4.67 
4.91 
4.51 
3.44 
3.50 


10  KPa  ksi 
32.2  6.36 


Compressive 
Strength 
ksi  KPa 


43.8  5  Mar  81 

40.1  20  Mar  81 

39.1  8  May  81 

23.8  25  Feb  81 

23.4  17  Jul  81 


E  6 

static  (xiQ 


Initial  Concrete  Elastic  Modulus 

/•~,«6  _  . ..  E ,  .  _ _ 6 


Date 

Tested 

6  May  82 

6  May  82 

6  May  82 

6  May  82 

6  May  82 


c  (X10^  psi) 

%  Increase 

5.59 

20 

5.89 

18 

5.43 

20 

4.54 

32 

4.54 

29 

Tota  1 

1'9 

119/5  =  23%  average  increase 

E  F. 

dynamic/ 'static  =  1.23 


Note:  All  cylinders  were  removed  from  fog  room  on  23  Feb  83. 

Rotating  specimens  produced  some  scatter  in  the  fnquency  and  damping 
values.  Some  of  the  specimens  were  slightlv  ovate. 


Table  4 


Frequency  and  Damping  Values 
for  First  10  Mode  Shapes 


4 


l 

! 

I 


< 


APPENDIX  A:  TRANSFER  AND  COHERENCE  FUNCTIONS 


and  coherence  function:  input-monol ith  7(1-7),  output-monolith  1(0-1);  el  495 


Transfer  and  coherence  function:  1-16,  0-9;  el  495 


14.0  16.0  18.0 


Transfer  and  co 


Transfer  and  coherence  functio 


Transfer  and  coherence  function:  1-22 


Transfer  and  coherence  function: 


Transf e 


and  coherence  function 


Transfer  and  coherence  function: 


is 


* 


ra 


APPENDIX  B:  EXTRACTING  MODAL 
PARAMETERS  USING  THE  HP5423A 
(Hewlett-Packard  1981  -  See  Measurement  Overview) 


1.  Once  the  data  stored  on  the  32  track  analog  system  are  digitized  by 
the  HP5423A  the  analyzer's  modal  analysis  capabilities  are  used.  From  the 
data,  transfer  functions  were  computed  representing  the  ratio  of  the  measured 
output  acceleration  to  the  measured  input  force.  A  total  of  70  transfer  func¬ 
tions  were  computed  from  the  available  data.  The  transfer  functions  (complex 
functions,  i.e.,  composed  of  both  real  and  imaginary  components)  are  used  to 
identify  the  resonant  frequencies.  The  resonant  frequencies  are  distinguished 
by  peaks  in  the  transfer  functions.  Once  the  resonant  frequencies  are  deter¬ 
mined  data  points  about  the  peak  are  curve  fitted  using  a  least  squares 
criterion.  Once  parameters  of  the  curve  fitting  function  are  determined  an 
estimate  of  the  damping  ratio  is  possible.  A  single  degree  of  freedom  (SDOF) 
function  is  used  to  fit  the  data  for  the  damping  ratio. 

2.  For  transfer  measurements  of  acceleration  to  force  peaks  occurring 
in  the  imaginary  part  and  zero  crossings  in  the  real  part  of  the  transfer 
function  mark  resonant  frequencies.  Also  the  phase  angle  between  the  real 
and  imaginary  parts  equals  ±90  degrees  at  a  resonant  frequency. 

3.  The  functions  shown  in  Figure  B-l  will  be  used  to  illustrate  how  the 
analyzer  is  used  to  compute  resonant  frequencies,  damping  ratios,  and  mode 
shapes.  The  transfer  function  from  the  response  point  on  monolith  13  is  used 
to  identify  the  resonance  at  5.94  Hz.  Looking  at  Figure  B-la,  a  peak  is 
observed  between  5.70  and  6.38  Hz  in  the  magnitude  of  transfer  function-- 
indicating  a  possible  resonant  frequency.  Next  any  of  the  three  remaining 
functions  can  be  studied.  Secondly,  the  phase  function  (phase  angle  of  the 
transfer  function  versus  frequency)  is  -90  degrees  at  about  5.9  Hz.  Thirdly, 
the  real  part  of  the  transfer  function  has  a  zero  crossing  at  about  5.9  Hz. 
And,  finally  in  the  imaginary  part  of  the  transfer  function  a  peak  is  observed 
at  5.94  Hz.  Therefore,  from  studying  the  four  functions  using  each  as  a  cri¬ 
terion  a  resonant  frequency  is  found.  If  other  transfer  functions  are  ob¬ 
served  peaks  are  found  at  5.94  Hz  also.  The  transfer  function  shown  was 


chosen  because  it  represents  a  point  on  the  dam  crest  with  the  largest 
response  recorded  at  5.94  Hz. 

4.  Once  the  resonant  frequency  is  found  11  data  points  in  the  transfer 
function  are  curve  fitted  using  a  least  squares  criterion.  The  data  used  in 
the  curve  fit  include  5  points  to  each  side  of  and  the  resonant  frequency 
points.  The  equation  is  in  the  following  form: 


H(S) 


r.  (S  +  2/ry  )  +  r  2rtf 
i _ so  r _ o _ 

S2  +  47tyoS  +  (27tyo)2  +  (27TfQ)2 


where  the  unknown  parameters  f  ,  y  ,  r.  ,  and,  r  are  determined  from 
the  least  squares  criterion  and: 

f  =  damped  natural  frequency  (Hz) 

Yq  =  damping  coefficient  (Hz) 
r^  =  imaginary  part  of  the  residue 
r^  =  real  part  of  the  residue 

H(S)  =  transfer  function  evaluated  at  S  =  j2nf  ,  where  f  is  frequency 
in  Hz 

j  =  V-i  _ 

12  2 

The  damping  ratio,  y^  ,  is  determined  by  Yq  ~  Yq/ v Y0  +  •  For 

instance  at  fQ  =  5.937  Hz  ,  yQ  =  0.2534  Hz  ,  therefore  Yq  =  (0.2534)/ 

V (0 • 2534 ) 2  +  (5.937)2  =  0.04262  or  4.262  percent  of  critical  damping  as 
shown  in  Table  4. 

5.  When  the  natural  frequency  and  the  damping  ratio  are  computed  resi¬ 
due  values  from  the  other  transfer  functions  are  determined.  The  residue 
values  are  proportional  to  the  magnitude  of  the  modal  shape  vectors.  For  ac¬ 
celeration  measurements  the  imaginary,  or  quadrature,  part  of  the  transfer 
function  is  used  to  compute  the  residues  for  the  other  response  points  as: 


r  =  r.  +  rr  -  2j  (2nY(j)  j  IJlKs)  ] ) 

(remember  r  is  a  complex  function)  where 

1  ( H ( S ) J  =  imaginary  part  of  the  transfer  function  evaluated  at 
m  S  =  j2nf 

o 

(Hence,  the  ruirne  "quadrature  peak  picking  method.") 


6.  The  residues  and  the  mode  shape  vectors  are  related  by  the 
expression : 

rk  =  <V  <V 

where 

rk  r  [rn(k)  r12^k)  •  •  •  1 

[r21(k)] 

[*] 

[*] 

[•) 

(where  r^(k)  is  the  residue  at  response  point  1  due  to  input  at 

response  point  2) 

U,  =  n  dimensional  vector 
k 

=  scaling  constant 

t  =  denotes  transpose  (refer  to  matrix  operations) 

7.  Hence,  when  all  the  residues  are  calculated  the  analyzer  uses  a 
routine  to  compute  and  scale  mode  shape  vectors  for  each  response  point.  The 
mode  shape  vectors  are  combined  and  define  the  mode  shape  for  the  portion  of 
the  structure  marked  by  the  response  points.  The  views  of  the  mode  shapes  in 
Figures  13  to  22  are  the  resulting  shapes  from  the  combined  mode  shape  vectors 
(one  vector  per  point  having  both  magnitude  and  direction). 
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